Observation of a Non-local Optical Vortex 
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We report the observation of an optical vortex in the correlations of photons produced from 
spontaneous parametric down- conversion. The singularity appears in a non-local coordinate plane 
consisting of one degree of freedom of each photon. 
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The spatial correlations of photons have received much 
attention recently due to the interest in quantum infor- 
mation and fundamental aspects of quantum mechanics 
E, i, i, i, S, . An important issue in this field is the in- 
vestigation of spatial entanglement through correlations 
of the orbital ang ular momentum (0AM) of photons 
0, 0, S, E El, [ia. In optics, 0AM arises from a phase 
singularity in the transverse profile of a beam of light 
around which the phase fronts twist, and for this reason 
these light fields are also known as optical vortices or 
screw dislocations [l3|. The motivation for the study of 
optical vortices ranges from fundamental interest to prac- 
tical applications. For example, 0AM of optical vortex 
beams has been used to rotate micro-particles flij . 
and to produce coherent atomic vortices in Bose-Einstein 
condensates [l6| and in cold ato ms [TtI . The 0AM of sin- 
gle photons has been measured [l8| and has been pro- 
posed as a carrier of high-dimensional quan tum infor- 
mation in quantum cryptograj 3hy [H, 0, Hll and quan- 
tum communication 



23l . |2j]. Furthermore, coher- 



ent transfer of 0AM between light and matter opens the 
possibility of storage and retrieval of high-dimensional 
quantum information [25|. 

Here we report the experimental observation of a non- 
local optical vortex associated with a phase singularity 
which manifests in the spatial correlations of two entan- 
gled photons. The vortex appears in a two-dimensional 
plane of the two-photon phase space that is composed of 
one spatial degree of freedom of each photon. The non- 
local character of this phase singularity makes it distinct 
from the usual 0AM entanglement of pairs of photons, 
in which the pump laser or down-converted photons are 
prepared or projected onto a light mode possessing 0AM. 

A common example of beams with optical vortices are 
the Laguerre-Gaussian (LG) modes, which are a set of 
field modes with well-defined angular momentum. The 
LG modes are parametrized by the radial and azimuthal 
quantum numbers p and m, respectively. A photon de- 
scribed by an LG mode carries an 0AM equal to mh. The 
LG beams are usually created with holographic masks 
[2^ or mode converters 27], the latter of which exploit 
the fact that a LG mode can be described as a superpo- 
sition of Hermite-Gaussian (HG) modes. For example, a 
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FIG. 1: (Color online) a) Intensity profile of a first-order 
Laguerre-Gaussian mode, b) Young's interference pattern 
from an infinitesimal double slit, represented by the vertical 
lines in a). The fringe patterns is shifted, due to the phase 
structure of the Laguerre-Gaussian beam. Along the horizon- 
tal dashed lines, the relative shift is half a period. 



first-order {jp = 0, m = ±1) LG mode can be written as 
a linear combination of HG modes [27 1 



{x,y) = A[ui{x)uo{y) ± iuo{x)ui{y)] , (1) 



where the HG modes lA^f^{x^y) = Uk{x)ui{y) have 
been decomposed into the one-dimensional HG functions 
Un{x)^ given by the field amplitude 
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Here Cn is a normalization coefficient, Hn{x) is the n*^- 
order Hermite polynomial and the parameters R{z)^ w{z) 
and (p{z) are the radius of curvature, beam waist, and 
Gouy phase, respectively. We note that uo{x) is the usual 
one-dimensional Gaussian function. 

Let us now consider transverse spatial degrees of free- 
dom of an entangled pair of photons. For simplicity, we 
use dimensionless position p = (px^Py) and wave vector 
q = {Qx^Qv) variables. In position space, the two pho- 
ton wave function produced by spontaneous parametric 
down-conversion (SPDC) can be written as [10, 28] 



^(pi,p2) = ^(pi + P2)r(pi -P2), 



(3) 
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FIG. 2: (Color online) Experimental setup. A non-local opti- 
cal vortex is produced in the correlations of entangled photons 
produced by SPDC in a BBO crystal. The two-photon state is 
engineered using the cylindrical lens to focus the pump beam 
and the Gaussian transmission masks (G) to spatially filter 
the down-converted photons. 

where £{p) is the field amplitude of the pump beam and 
r(p) is the phase matching function. Let us approximate 
the phase matching function by a Gaussian of width 5: 
r(p) = Uoo{p^S) = uo{px^S)uo{py^S)^ and consider that 
the pump laser is given by a first-order HG mode with 
width cr: £{p) = Uio{p,cr) = Ui{pj,, cr)uo{py, a), where 
notation has been modified to include the widths explic- 
itly. Then, the two photon wave function is given by 

(Pi, P2) = ^io(pi + P2, cr)Uoo{pi - P2, S). (4) 

Choosing a = and using the fact that the HG func- 
tions are eigenstates of the Fourier transform, the two- 
photon state can also be described in a mixed position- 
momentum representation by the wave function 

^ (Pi, <72) =C [ui{pxi)uo{qx2) + iuo{pxi)ui{qx2)] 

X uo{pyi)uo{qy2), (5) 

where C is a normalization constant. The x-component 
of the wave function is equivalent to the LG mode U^q 
given in Eq. ([T]), where pxi and qx2 play the role of 
the X and y variables. Thus, the two-photon field in the 
Pxi X qx2 plane presents an optical vortex. The prob- 
ability distribution P{pxi^qx2)^ obtained by integrating 
IV^ {PI1Q2) ? over pyi and qy2, is 

P{pxi,qx2) = \uo{pxi)ui{qx2) ^ iui{pxi)uo{qx2)\'^ , 

(6) 

which corresponds to the modulus squared of the LG 
mode given in Eq. ([1]), defined in non-local coordinates 
Pxi and qx2' 

Using photons produced from spontaneous paramet- 
ric down conversion (SPDC), we generated a two-photon 
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FIG. 3: (Color online) a) Coincidence count distributions cor- 
responding to the probability distribution P{pxi,px2) are de- 
scribed by a Her mite- Gaussian mode, b) Coincidence count 
distributions corresponding to the probability distribution 
P(pxi,qx2) are described by a Laguerre- Gaussian mode. The 
zero intensity region at the center is due to the undefined az- 
imuthal phase at the origin. The arrows show the positions 
of detector 1 used to obtain FIG.Hl 



state accurately described by Eq. ([5]) and observed a 
non-local optical vortex. The experiment is shown in 
FIG. [21 A 441.6 nm c.w. HeCd laser is used to pump 
a 1cm long BBO crystal cut for type II phase match- 
ing. The degenerate entangled photon pairs were sent 
to single photon detectors outfitted with interference fil- 
ters (40nm FWHM) centered at 884 nm. Generation of 
the state ([5]) consists of three steps. In the first step, 
a 125/im thick microscope slide is inserted into half of 
the laser beam. The angle of the slide is adjusted so 
that there is a tt phase shift between the two halves of 
the Gaussian beam. After spatial filtering, the laser pro- 
file is approximately described by a HG mode U^q^ . In 
the second step, the pump laser is focused at the cen- 
ter of the BBO crystal using a cylindrical lens with focal 
length 33mm. This guarantees that the widths a and S 
of the pump beam profile and phase matching function 
F are approximately the same. In the third step, the 
wave vector distributions of the down-converted photons 
at the center of the crystal are mapped onto intermediate 
planes through optical Fourier transforms (/i = 100mm). 
The field profiles at the intermediate planes are spatially 
filtered using transmission masks G with a Gaussian den- 
sity profile, which guarantees the Gaussian shape of the 
function F. The two-photon state at the intermediate 
planes just after the masks G is given by Eq. ([5]), where 
pi and p2 are the transverse position coordinates. 

To verify the presence of the HG mode in the entan- 
gled two-photon wave function, the intermediate planes 
are mapped onto the detection planes using an imag- 
ing system composed of two confocal lenses with focal 
lengths /2 = 150mm and fs = 50mm. The Pxi x px2 cut 
of the phase space distribution is measured by scanning 
the two detectors in the vertical direction, resulting in 
a two-dimensional (2D) array of coincidence counts that 
is proportional to the detection probability P{pxi^ Px2), 
as shown in FIG. [3] a). In this figure one can see 
that the coincidence count distribution along the diag- 
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FIG. 4: (Color online) Interference fringes when the pump 
beam is a) a Hermit e- Gaussian beam and b) a Gaussian beam. 
In both figures the sets of fringes are obtained by fixing detec- 
tor 1 at positions pxi = 0.2 (red circles) and pxi = 0.4 (blue 
squares) and scanning detector 2. Solid lines correspond to 
curve fits of Eq. (0. 



onal, which displays correlations coming from the func- 
tion E{pxi + Px2)^ is indeed given by first-order Hermite- 
Gaussian mode. The anti-diagonal direction corresponds 
to correlations coming from the function T{pxi — ^^2), 
and is approximately described by a Gaussian function. 
Next, the Fourier transform of the intermediate plane G 
of photon 2 is performed using a lens with focal length 
/4 = 250mm. Since the detection plane lies in the focal 
plane of this lens, the field distribution in the detection 
plane corresponds to the Fourier transform of the plane 
G up to a quadratic phase term, which does not con- 
tribute to the intensity distribution. FIG. [3] b) shows 
the 2D array of coincidence counts corresponding to the 
probability distribution P{pxi^ Qx2)^ which displays the 
same doughnut shape that is characteristic of the inten- 
sity profile of the LG field mode, as shown in FIG. [1] 
a). The coincidence distribution is not a perfect repli- 
cation of FIG. [1] a), due to noise which is most likely 
due to slight misalignment of the lenses, the Poissonian 
photon count statistics, and the finite width of the detec- 
tor apertures (20/im for p measurements and 50pm for q 
measurements). Despite this noise, the doughnut shape, 
in particular the zero intensity region at the center, is 
clearly visible. 

Even though the coincidence profile in FIG. [3] b) re- 
sembles the intensity profile of a LG mode, in order to 
verify the existence of a phase singularity it is necessary 
to characterize the phase dependence. A simple method 



to identify the azimuthal phase dependence of an opti- 
cal vortex is through interference in a typical Young's 
double slit experiment [29], as illustrated in Fig. [T]b). 
The intensity of the Fraunhoffer diffraction pattern of a 
Laguerre-Gaussian beam incident on a double slit, ori- 
ented to produce fringes in the x direction, is given by 



I{x,y) = R{x,y) 



1 



a{x) + A(l){y) 



(7) 



where R{x^y) is a Gaussian function, a{x) is the usual 
phase difference due to propagation from the each slit to 
the detection plane and A(j){y) = 0i(^) — (/>2(^) is the 
phase difference of the input LG field at each slit. We 
emphasize that the interference fringes now depend on 
both X and y directions. The azimuthal phase structure 
results in an oblique set of fringes as illustrated in FIG. 
[T]b), in contrast to the usual vertical linear fringes which 
appear in Young's interference of a plane wave or Gaus- 
sian mode. Therefore, the fringe pattern obtained by 
scanning in the x direction for different y positions will 
be phase-shifted from one another. Choosing the y po- 
sitions represented by the dashed lines in FIG.[T]b), the 
shift can be close to half a period. 

To characterize the azimuthal phase structure, photon 
1 was measured using the imaging system, while a double 
slit (200/im slits, 100 pm separation) was placed in the 
focal plane of lens /2 of photon 2. The field distribution 
in this plane corresponds to the wave vector distribution 
of the field at the plane of the mask G. Lens fs was 
then used to map the Fourier transform of the field dis- 
tribution along the double slit onto the detection plane 
of detector 2, corresponding to the Fraunhofer diffraction 
regime. Detector Di was placed in two positions (0.2mm 
and 0.4mm) indicated by the arrows in FIG.[3]b), while 
detector 2 was scanned linearly in the transverse plane, 
corresponding to the horizontal axis in FIG. [3]b). The 
result was the set of interference fringes shown in FIG. 
m a) . A relative shift of approximately half a period be- 
tween the two sets of fringes is clearly observed. For 
comparison, we removed the glass microscope slide from 
the pump laser, obtaining a Gaussian profile beam. The 
corresponding interference fringes are shown in FIG.[4]b). 
In this case no shift of fringes is observed, demonstrating 
that the features in FIG. [H a) are indeed a consequence 
of the azimuthal phase dependence of the spatial correla- 
tions. The asymmetry of the fringes in FIG.[4]b) is due 
to a small distortion in the shape of the Gaussian beam 
in the non-local plane. 

It is interesting that the non-local vortex appears in the 
coordinates of entangled photons, which suggests that 
there may be a relation between entanglement and iso- 
lated phase singularities. We will now show that appear- 
ance of an isolated phase singularity implies that a bipar- 
tite pure state is entangled. Consider the wave function 
(j){x^y) = \R{x^y)\ex.-p[ih{x^y)]^ where x and y are co- 
ordinates for two modes. Let us choose the vortex such 
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that h{x^y) is undefined at the origin. For (j) to be separa- 
ble requires h{x,y) = hi{x) + h2{y), which imphes that 
hi and/or /12 is undefined at x = and/or = 0, re- 
spectively. But this implies that y) is singular at all 
points along the line x = and/or ^ = 0, which contra- 
dicts the original assumption of an isolated singularity. 
Thus, in the present case the observation of the non-local 
vortex implies in quantum entanglement. 

Let us briefly discuss our results in terms of propa- 
gation of the two-photon state, considering one spatial 
dimension of each photon. At the source the spatial pro- 
file of coincidence counts is given by the diagonal HG 
mode shown in FIG.[3ti). After photon 1 passes through 
an imaging system and photon 2 through a Fourier trans- 
form system, the coincidence profile is given by the LG 
mode shown FIG. [8)3). This propagation is analogous to 
the mode conversion of a classical HG beam to a LG 
beam, achieved through introduction of an astigmatic 
propagation region [27[. Here the astigmatism occurs 
between the propagation of photons 1 and 2, and cor- 
responds to the photons propagation through different 
lens systems: the imaging system of photon 1 and the 
Fourier system of photon 2. In this sense, propagation 
of the two-photon state from the Gaussian masks to the 
detection plane can be considered as a non-local mode 
conversion. 

In conclusion, we have observed a phase singularity in 
a non-local coordinate plane of a pair of entangled pho- 
tons. The phase singularity appears as a non-local optical 
vortex in the spatial correlations of the entangled pho- 
tons. Though we have shown the existence of a first-order 
vortex, it should be possible to create higher-order vor- 
tices, or more complicated structures. We also note that 
these non-local singularities should appear in any pair 
of entangled harmonic oscillators, such as modes of the 
electromagnetic field or the vibrational modes of trapped 
ions. One might envision using these types of vortices to 
encode quantum information non-locally in two photons 
or two atomic clouds, or to produce a two-photon angular 
momentum transfer effect in two-photon absorption. 
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